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ABSTRACT 

We study the effect of a white dwarf on the spin-down of a cataclysmic variable system using a 
three-dimensional magnetohydrodynamic numerical model. The model includes the stellar corona, 
the stellar wind, and the WD mass and magnetic field. The existence of the WD modifies the system 
spin-down by physically blocking the stellar wind, restructuring the wind, channeling the wind towards 
the WD surface, and by modifying the shape and size of the Alfven surface. The combination of 
these processes differs among a set of simple test cases, and the resulting angular momentum loss 
rates vary by factors of 2-3, and by factors of two relative to a test model with a single M dwarf. 
While the model employs some simplifications, the results suggest angular momentum loss schemes 
currently employed in cataclysmic variable studies do not require drastic revision. Insights are also 
gained on wind accretion. We find that efficient accretion switches on quite rapidly with decreasing 
orbital separation. Accretion rates depend on magnetic field alignment and should be modulated 
by magnetic cycles on the M dwarf. For particular values of white dwarf magnetic field strength, 
an efficient syphoning of coronal plasma from the inward facing M dwarf hemisphere occurs. Wind 
accretion rates are expected to vary by factors of 10 or more between fairly similar close binaries, 
depending on magnetic field strengths and orbital separation. 
Subject headings: stars: magnetic field — stars: coronae — binaries: close 



1. INTRODUCTION 

The orbital and mass transfer evolution of cataclysmic 
variables (CVs) and their pre-CV immediate antecedents 
is though to be largely governed by angular momentum 
loss through gravitational radiation and to the magne- 
tized wind of the unevolved M dwarf component ( |Kraft| 
et al. 1962 Verbunt & Zwaan 1981 ). Angular momentum 
loss through gravitational radiation is well-understood 
for the stars in CVs, but there is still no comprehensive 
theory of spin-down through magnetized winds that is 
expected to dominate by an order of magnitude or more 
for periods of > 3hr, above the CV period gap flRap- 
paport et al. 1983 Andronov et al. 2003 Ivanova et al. 
20030 . 

Magnetic braking depends on the wind mass loss rate 
(MLR) and the large-scale stellar magnetic fiel d (e.g. We- 
ber fc Davisjl967||Mestel|1968| pCawaler|1988| ). Both are 
extremely difficult to measure for late- type stars and are 
especially uncertain at the very rapid rotation rates of 
close binaries where magnetic proxies such as X-ray emis- 
sion show saturation and supersat uration effects that ap- 
pear to affect spin-down rates (e.g.|Andronov et al.|2003| 
Ivanova & Taam|2003| |Barnes|2003[ I Wright et al.||2011). 

^ I f | I f | o II V 

Moreover, it is generally assumed that magnetic braking 
follows the same prescription as f or single stars and that 
the WD makes no difference (e.g. |Rappaport et al.|1983 



Andronov et al. 2003) [ivanova fc Taam|2003| |Davis et al 
2008). While this might be reasonable for WDs with 
weak magnetic fields that have no direct influence on the 
magnetospheric structure of the system, it is not clear 
that such a situation applies to all systems. In particular, 
WDs can have surface equatorial magnetic fields ranging 
up to 10 s G that would be expected to have some con- 
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siderable influence on the magnetically driven wind of a 
close late-type companio n (Schmidt et al. 1986; Ha meury 
et al.|1987||Li et al.|1994b|a[|1995||Li fc Wickramasinghe 



1998) 



Semi-empirical magnetohydrodynamic (MHD) models 
of the solar wind have now reached a level of realism 
such that they are able to repr oduce quite accura tely 
salient solar wind properties (e.g. Cohen et al.||2007| ). In 
previous work, we applied such models to other stars, in- 
cluding rapid rotators and st ars harbouring giant planets 



(Cohen et al. 2009 2010b|a ). Here, we apply the same 



modelling techniques to begin to tackle the problem of 
stellar winds and angular momentum loss in close bina- 
ries and examine the influence of the companion WD on 
magnetically-driven stellar winds. In this pilot study, we 
consider a synchronized system, such that the WD and 
late-type star are locked in phase and rotating with the 
orbital period. 

2. SIMULATION 



This work follows from Cohen et al. (2010c), who stud 



ied the effect of a close-in planet on the stellar angular 
momentum loss to the stellar wind as a function of the 
planetary semi-major axis. They found that the stel- 
lar angular momentum loss rate (AMLR) decreases once 
the planetary Alfven surface, which is the "obstacle" the 
magnetized corona and wind feel, starts to interact with 
the stellar Alfven surface. Here we focus on two other im- 
portant factors that should affect the AMLR of a close 
binary system: the mass and magnetic field of a WD, 
which are much greater than those of the planet. 

We use the solar corona numerical model described by 
Cohen et al. (20071 120081), which is based on the generic 
BAl'S -R- U S MTlTJ mocTel ( |Powell et al.|1999||Tl3tJllFaT 
2012|. The model solves the set of conservation laws 
for mass, momentum, magnetic induction, and energy 
with some numerical resistivity, and is driven by surface 
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magnetic field maps that are used both to determine the 
initial (potential) magnetic field distribution as well as 
to scale the boundary conditio ns on the stellar ph oto- 
sphere. We refer th e reader to |Cohen et al.| (|2009| and 



Cohen et al. ( 2010c[ ) for further details and the methods 
we employ to calculate the model MLR and AMLR. 

For our CV-like system, the donor star is at the center 
of the simulation box and the WD at f w d = (3i?*,0,0). 
The numerical MHD code provides a self-consistent wind 
solution for the star, in which the WD is represented by 
an additional boundary condition for the density, tem- 
perature, and magnetic field. In addition, the gravi- 
tational potential due to the WD mass is included in 
the MHD equations. The very small scale of the WD 
and very high surface magnetic field can cause signifi- 
cant numerical difficulties for the MHD model. How- 
ever, the effect of the WD on the system, even out to 
relatively large radial separations, depends only on its 
central coordinates, mass and dipole moment, and not 
on the details of its body itself. We therefore truncate 
the simulation domain so that the WD boundary (i.e., 
the WD size) is 0.2i?* instead of roughly an Earth size 
body (~ O.Oli?*). This enables us to have high resolu- 
tion (small grid size, Ace) around the WD and lowers the 
magnetic field strength, B, near the boundary, which is 
a key constraint on the simulation time-step. 

We adopt a stellar radius and mass of = 0.5Rq and 
M* = 0.5M Q , for the donor (the subscript denoting so- 
lar values) and Mwd = O.75M for the WD. The star is 
assumed spherical: Roche lobe geometry and associated 
distortion of the secondary star is neglected. The donor 
and WD magnetic fields are dipolar with equatorial field 
strengths of f?* = 5 G and Bwd — 777 G, respectively. 
The latter is equivalent to an equatorial field strength of 
B w 10 6 G at the actual WD surface. Due to numerical 
limitations, both the stellar and the WD magnetic fields 
are an order of magnit ude lower than the fields recently 
obser ved in M dwarfs (Phan-Bao et al. 2009 Morin et al 



fects of WD gravity, field strength and alignment, and 
orbital separation. The parameters are summarized in 
Table [I] We also computed a reference model for the 
M dwarf alone to judge the influence of the WD on the 
wind solutions compared to the single star case. 

The simulation box extends from — 45-R* to +45i?* in 
each dimension, and a non-uniform grid is constructed so 
that high-resolution (Ax « O.Oli?*) is obtained around 
the two bodies and near the equatorial plane. We com- 
puted the resulting MLR and AMLR normalized to the 
values for the reference single star case for which we ob- 
tained a total MLR M = 4.1 • 10 12 [g s _1 ] and an AMLR 
j = 7.4- 10 30 [g cm- 2 s~ 2 }. 

3. RESULTS & DISCUSSION 

3.1. Model Insights 

In our simulation, the coronal magnetic field and the 
stellar wind are computed self-consistently under the dis- 
ruptive influence of the WD. The steady-state wind solu- 
tion for each set of parameters defines the shape and size 
of the Alfven surface (the collection of points where the 
wind becomes super- Alfvenic). The stellar MLR is sim- 
ply the integral of the mass flux over this surface. The 
AMLR is the integral of the product of the mass flux at 
the surface (that applies the spin-down torque) with the 
local radius squ ared times the rota tion rate (see equa- 



tions 5 and 6 in Cohen et al. (2009)) 



As we shall see, the WD magnetic field results in sig- 
nificant opening, closing and distortion of the M dwarf 
field down to the stellar surface in ways that differ ac- 
cording to dipole alignments, WD field and orbital sep- 
aration. The MHD wind model employs the empiri- 
cal Wang-Sheeley-Arge (WSA) relation th at relates wind 
speed to the magnetic field ex pansion rate (|Wang fc Shee 



Iey|p90 



of the fie 



Arge & Pizzo 



d, the slower 



2000): the greater the expansion 



2010 1 , an d those used m previous stud ies of magne tic C V 
syste ms ( |Li et al!||1994b|a| |1995| Li & Wickramasinghe 
1998]) . It is not yet clear ho w stellar winds scale wit 
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he resulting wind. The different 
magnetic field topologies in the cases studied are then 
naturally expected to affect the resulting wind proper- 
ties. 

Figure [T] shows the density structure and the magnetic 
field topology in the region between the donor and the 
WD for all six cases. There is a region of high density 
near the star, as well as a ring around the WD, where 
channelled plasma is trapped within closed field lines. In 
cases A and C, the WD and donor dipoles are aligned, 
and a pressure balance point is seen where the two fields 
"collide". In the anti-aligned cases (B, D, E, and F), a 
magnetic (neutral) X-point is seen where the stellar and 
WD fields reconnect and many field lines connect the 
two bodies. This leads to an accretion of stellar wind 
plasma onto the WD and an expansion of the origin of 
the high-density wind towards the donor poles. 

The effect of the WD mass is clearly seen when compar- 
ing cases A and C, and B with D. The WD mass focuses 
the stellar wind to accrete along field lines, increasing the 
density in the region between the objects. The pressure 
balance point (cases A and C) and the magnetic X-point 
(cases B and D) both move towards the WD when its 
mass is included. 

3.2. Mass loss, angular momentum loss and accretion 

Table [l] summarizes the MLRs, M, and AMLRs, J, for 
We performed six case studies, A-F, probing the ef- the different test cases normalized to the respective single 



e stellar field strength (e.g. Cohen et al. 2009, 2 010a[ ), 
but the magnetic balance of our system — the ratio of M 
dwarf and WD field strengths — approximates a magnetic 
CV with a WD field B WD « 10 7 G. 

CVs typically comprise a tidally locked, synchronously 
rotating secondary, and a non-synchronously rotating 
WD surrounded by an accretion disk. We treat a sim- 
plified synchronous, diskless case here as a first step; 
treatment of the more complex asynchronous system re- 
quires a much more computationally demanding time- 
dependent solution beyond the scope of the work in 
hand. We omit centrifugal effects and investigate the 
non-rotating, tidally-locked scenario, neglecting any rel- 
ative azimuthal motion between the WD and the coronal 
plasma. This synchronously rotating WD and secondary 
is more analogous to the situation in AM Hcr-types (de- 
spite our adopted WD magnetic field strength, which is 
an order of magnit ude below th at r equired to synchronize 
the W D rotation; |Joss et al.| e.g. [19791 |King fc White^ 
hurst e.g. 1991 Wickramasinghe fc Wu| e.g^ 1991D and in 



magn etic pre-OVs, or "pr e-polars" (e.g. Schwope et al 
2002] |Schmidt et al.|[2005 | 
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M star values, Mq and Jo. The last two columns list the 
mass accretion rate, M acc , onto to the WD (normalized 
to Mo) and the total unsigned radial magnetic field at 
the Alfven surface, B Q (the total magnetic flux opened 
by the stellar wind). This flux is not normalized to the 
reference case since it is dominated by the strong WD 
field. 

The distribution of the mass flux and the shape of the 
Alfven surface for the six cases are shown in Figure [21 
The reference case without the WD is similarly displayed 
in Figure [3j In general, the MLR is governed more by the 
plasma density than its speed. Wind velocities computed 
from the model tend not to vary by more than a factor 
of 2 or so, as is the case for the solar wind, whereas the 
wind density is much more inhomogeneous. The most 
notable feature in all cases is the reduction of the MLR 
in the hemisphere occupied by the WD. 

There is a 70-85% reduction in MLR for cases A and 
B for a massless WD (equivalent to a planet with a very 
strong field!). The MLR for aligned dipoles in case A 
is approximately twice that for the anti-aligned case B. 
For anti-aligned dipoles there is a large region of closed 
field between the stars that prevents wind escape. Field 
lines pulled toward the WD increase the field expansion 
factor at the poles, weakening the wind. For the aligned 
case A, the closed field region is much smaller, and field 
lines are compressed with smaller expansion factors near 
the poles. The WD accretion rate for these unphysical 
massless WD cases is actually negative and is caused by 
the boundary conditions serving as an evaporating mass 
source. This boundary effect is removed when the WD 
mass is turned on. 

The reduction in MLR is not as extreme — about 50% 
for both cases C and D — when the WD mass is included, 
though for case C there is a large increase in mass ac- 
creted by the WD. When moving the WD further out 
(case E), the reduction of MLR is only 10%. Here the X- 
point is further from the donor and its closed field corona 
covers a larger area, shutting off some of the wind. Cou- 
pled with a reduced influence of the WD gravity close 
to the donor surface, the larger orbital separation leads 
to a cut in accretion by an order of magnitude. In the 
weaker WD field (case F), the MLR is reduced by 70% 
due to the very efficient magnetic coupling that develops 
between the donor star and WD. The X-point for this 
case is much closer to the WD surface than for case D, 
leading to smaller field opening angles and efficient "sy- 
phoning" of coronal plasma from the donor. 

The AMLR is governed by the total MLR through the 
Alfven surface and by the radial extent of the surface. 
The Alfven surface is inflated compared with the single 
star case, with the inflation being greater for the mag- 
netically aligned fields in cases A and C. For these cases, 
the inflation of the surface is more significant than the 
decrease in MLR, so that the AMLR is 30% (case A) and 
60% (case C) higher than the reference case. In cases B 
and D, the inflation of the surface is less significant than 
the reduction in MLR, so that the overall AMLR is only 
about 10% lower than the reference case. It is interest- 
ing to note that by only moving the planet from 3i?* to 
4.5i?* in case E this trend is reversed and the AMLR is 
25% higher than the reference case. Despite the signifi- 
cant increase in Alfven surface size, the AMLR in case F 



is 50% lower than the reference case. This demonstrates 
the importance of the details of the Alfven surface on 
the AMLR: the extended protrusions in the surface in 
the left hemisphere where the wind density is reasonably 
high that are particularly prominent in case D are absent 
in case E. Despite the large equatorial extension of the 
surface in case F in the right hemisphere there is little 
effect on the AMLR because of the very low wind density 
in that region. 

4. CONCLUSIONS AND IMPLICATIONS FOR 
BINARY EVOLUTION 

We note that there are several ways in which the work 
presented here could be extended and improved: mod- 
eling of an asynchronous rotation of the star and WD; 
investigation of field alignments inclined relative to the 
system spin axis; inclusion of more complex M dwarf sur- 
face magnetic field distributions; treatment of magnetic 
field strengths stronger by at least an order of magni- 
tude to probe more rigorously the parameter space of 
magnetic CVs; inclusion of an accretion disk for lower 
field systems. 

With this caveat in mind, the salient point from the 
numerical models is that both the wind-driven MLR and 
AMLR from an M dwarf can be significantly affected by 
the presence of a WD companion — by factors of 2-3 for 
the limited parameter range studied. The alignment of 
the component magnetic fields can also be important, 
with differences between aligned and anti-aligned fields 
amounting to nearly a factor of two in AMLR. While 
the differences between the single star and CV-like mod- 
els are significant, they are not drastic. We conclude 
from this pilot study that gross changes to AMLR for 
spin-down driven by M dwarf winds appear not to be 
required, especially when noting th at different AMLR 
approaches adopted to date (see, e.g. Knigge et al.|2011 1 
already harbor uncertainties as large as, or larger than, 
the differences between our various models. 

The possibility that a strong magnetic field on the 
WD of a close binary might influence the magnetically 
driven wind of the la te-type companion was suggested by 
Schmidt et al. ( 1986[ ), who estimated that the 10 G fields 
of WD primaries m the magnetic AM Her-type systems 
could increase the magnetic lever arm and the AMLR 
through the secondary w i nd by fac t ors of 2- 5. Wick 
ramasmghe fc Wu| ( |lt)94l ; |Li et al.| ( |1994bla| |1995| ); |L 
& Wickramasinghe] ( |1998p calculated AMLR and othei 
properties of magnetic C V s by superimposing the poten- 
tial fields of the two components, and using supersonic 
flow and density profiles which depend on the gravita- 
tional potential of the two objects. They argued that 
strong WD magnetic fields could inhibit the donor wind 
outflow, reducing open field regions and instead captur- 
ing the secondary wind and reducing the magnetic brak- 
ing relative to single stars. They predicted that the wind 
would be completely captured for fields of several 10 7 G. 

The MHD wind model employed here represents a large 
improvement over the earlier methods that were based on 
a potential field approach. The dependency of the wind 
properties (velocity and density) on the field structure, 
as well as the self-consistent evolution of the field struc- 
ture as it is dragged by the magnetized wind, cannot be 
treated by potential fields. The drawback to our models 
is that we are not yet able to probe the 10 7 -10 8 G fields 
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investigated in the earlier studies. This would be valu- 
able for testing whether the mass loss rates from the M 
dwarf under such an external field influence are indeed 
similar to that for a single star. However, there are three 
aspects of our results of relevance for the wind accretion 
phase of systems with weaker WD fields. 

Firstly, the accretion rate depends quite strongly on 
magneti c field alignment, as pointed out for stronger 
fields by |Li &: Wickramasinghe| ( [l998[ ) , with a factor of 5 
increase from case C (aligned dipoles, low accretion rate) 
to case D (anti-aligned dipoles, higher rate). During a 
wind capture phase such as is thought to characterize 
the "pre-polars" — magnetic systems not yet filling their 



Roch e lobes (e.g. |Schwope et al. 2002 Schmidt et al 



2005|) — we might then expect to see accretion signatures 



vary with any magnetic cycles of the donor star. Such 
modulation could provide novel insights into magnetic 
cycles in stars near the fully-convective limit. 

Secondly, we see accretion is reduced by an order of 
magnitude when the orbital separation is increased from 
3-R* to 4.5i?* (case E vs. case D). At face value this 
suggests that wind accretion effectively "turns on" as 
the system evolves to shorter periods much more quickly 
than would be expected by purely hydrodynamic accre- 
tion. In the Bondi & Hoyle (1944) recipe, the accretion 
rate varies as the inverse square of the orbital separation 



(see, e.g., Matranga et al. 



2011) implying only a factor 



Thirdly, we find a large increase in mass accretion rate 
for case F with a weaker WD field of 10 5 G. For this 
anti-aligned case, effective magnetic coupling between 
the stars develops, and what would have been closed- 
field corona is effectively syphoned off onto the WD. The 
accretion rate is equivalent to the entire wind of the sin- 
gle star case. This "magnetic syphon" case then mimics 
the accretion rate of the total wind capture scenario of 
the pre-polars. The accretion rate is r eminiscent of the 
2 x 10 ~ 13 M(7) yr _1 recently found by Matranga et al. 
(2011) for QS Vir, a non-magnetic pre-CV that does 



of 2 difference in rates between cases D and E 



not appear to fill its Roche lobe. Such an accretion rate 
is otherwise difficult to understand under hydrodynamic 
wind capture, appearing an order of magnitude too high 
for likely M dwarf wind mass loss rates. 

The variation in wind accretion rates by factors of more 
than 10 within the fairly narrow range of model parame- 
ters investigated here suggests that wind accretion rates 
among pre-CVs could cover a wide range of values. 
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from NASA-ESS, NASA-ESTO-CT, NSF-KDI, and 
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TABLE 1 

Summary of simulated test cases 





Case 


a [R+] 


Bwd [G] 


M/Mo 


J/Jo 


M acc /M 


Bo [G] 


A 


- aligned B wl i, massless WD 


3 


10 6 


0.31 


1.39 


-0.02 


60.24 


B 


- flipped B w d, massless WD 


3 


10 6 


0.14 


0.90 


-0.03 


61.15 


C - 


aligned B w< i> with WD mass 


3 


10 6 


0.51 


1.61 


0.09 


46.43 


D - 


flipped B w( j t with WD mass 


3 


10 6 


0.48 


0.92 


0.55 


42.09 


E - 


flipped B w d, with WD mass 


4.5 


10 6 


0.89 


1.26 


0.05 


53.15 


F - 


flipped B wt i, with WD mass 


3 


10 5 


0.34 


0.48 


0.95 


23.33 



Case B 




Fig. 1. — The y=0 meridional cut zoomed to the region between the star (red sphere) and the WD (white sphere). Color contours are of 
number density and white streamlines represent magnetic field lines. The solid white line in the two left panels shows the location of the 
magnetic pressure balance of Case A, while it shows the location of the x-point for Case B in the two right panels. 




Fig. 2. — The y=0 meridional cut colored with mass flux distribution(in g s 1 cm 2 ). The Alfven surface is shown in solid black line, 
while the star and the WD are shown as red and white spheres, respectively. 
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Fig. 3. — The same display as in Figure[2]for the reference case without the WD. 



